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P I In this paper, we develop a non-Bayesian methodology to analyze the time- 

varying structure of international linkages and market efficiency in stock markets. 
■^ . We consider a non-Bayesian time-varying vector autoregressive (TV-VAR) model, 

and apply it to obtain the time-varying generalized impulse responses (TV-GIRs) 
f-H ' between the U.S. and Japanese stock markets. Our empirical results based on the 

C/j . TV-GIRs provide a new perspective that the stock market linkages and market effi- 

^H . ciency in the semi-strong sense change over time. We conclude that their behaviors 

^+H I correspond well to historical events of the international financial system. 
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^ : 1 Introduction 



There is much literature on international stock market linkage in 1990s and early 2000s. 
The literature on the topic can be categorized into two approach es. The first approac h 



employs a vector autoregressive (VAR) model (see, for example, lEun and ShimI (119891 ) 
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Jeon and Von Furstenberd (119901) . iMathur and Subrahmanvaml (Il990l ) , IChoudhrvl (119941 ) , 
and iTsutsui and Hiravamal (120041 ) ). These studies are concerned about the existence of 
international hnkage of stock prices, especially about the short-run relationships among 
stock price indices. The second approach employs a vector error correction (VEC) model 
and makes cointegration tests to exa mine a long-run e q uilibr iu m relationsh i p among the 
return s of stocks (see, f o r example. I Jeon and Chiand (Il991). IChan et al 



'4J. ll997h ). 



1992h. iKasa 



Lai et al.l (119931). ICorhav et al.l ( 11998h . Isiackman et al.l (11994h . IChung and Lhi 



(fl99i 

(119941 ). and IChoudhrvl (11994 119971 )). Some of them study both short-run and long-run 
relations by investigating cointegration first and then estimating a VEC model. They 
concerns a cointegration found to exist as a long-run equilibrium. Furthermore, some 
studies of the second group examine variance decomposion and impulse response functions. 
They have common interest in dynamic relationships of international stock market with 
ours in a sense. 

The preceding studies concerning linkage of stock markets employ a VAR model, 
conventional VAR one, structural VAR one or VEC one. They essentially suppose that 
parameters in the models are invariant over time. In other words, there supposed to be no 
structural change over the periods of their observation. This supposition is inappropriate 
when we consider that our contemporary financial system usually confronts with large 
or small perturbations, which sometimes cause very big structural changes such as the 
recent bankruptcy of Lehman Brothers. Whatever model is adopted, they lead to a flaw if 
one consider the coefficients of the VAR model as time-invariant within the observations. 
Such a model deals with the time- varying structure of international stock market linkage 
at most by taking subsamles or by adopting a rolling method. It cannot avoid arbitrary 
choice of the size of subsamplels and of the width of rolling window. Although several 
studies using VEC models and taking cointegration approach regard dynamic aspect of 
stock market linkages, they also fail to avoid the flaw. That is, their results are sensitive 
to the number of observations and the lag order of which VAR model is chosen. 

In this paper, we develop new methodologies to analyze time-varying structure of in- 
ter nationaHinkagesandjnarket efficiency in stock markets as practice for an univariate 
data in llto and Nodal (120121 ). We consider a non-Bayesian time- varying VAR (TV- VAR) 
model, and apply it to obtain time- varying generalized impulse response (TV-GIR) func- 
tion between the U.S. and Japanese stock markets. In contrast to the approaches of the 
preceding studies, our TV- VAR approach enable us to avoid their arbitrariness. We can 
shed new light on the international stock price linkages form the view point of dynamic 
change of international stock market. In particular, we focus on the degree of market effi- 
ciency and propagation of market fluctuations in the semi-strong sense by using TV-GIR 
function and long-run multipliers. Our empirical results provide a new perspective that 
the stock market linkages and market efficiency in the semi-strong sense change over time. 
We conclude that their behaviors correspond well to historical events of the international 
financial system. 

This paper is organized as follows. Section |2] presents the model and new methodolo- 
gies for our non-Bayesian TV- VAR model. Section |3] describes the data used. Section 
m shows the empirical results that the stock market linkage and market efficiency in the 
semi-strong sense change over time and that their behaviors correspond well to historical 
events of the international financial system. Section|5]is allocated to concluding remarks. 



2 Model and Methodology 

This section shows our model and methodologies to elucidate the time- varying structure 
of international linkage and market efficiency of stock markets. Some preceding studies 
concerning the subject have employed VAR models portioning whole sample into several 
subsamples or rolling subsamples with fixed window width. In place of such an approach, 
we develop a new method. We first examine whether the VAR estimates change over time 
or not by using Hansen's (1992) parameter constancy test. When the null of Hansen's 
(1992) parameter constancy test is rejected, we employ a non-Bayesian TV- VAR model 
as the alternative of Hansen's test assuming that the coefficients follow a random walk. 
Then, we estimate our TV- VAR model as a state space model by regarding it as a huge 
linear regression model. Considering the state space model such as patchwork of locally 
stationary models, we can analyze the the dynamic stock market linkage for period by 
period as well as market efficieny. In particular, the TV- VAR coefficients estimated 
allow us to obtain TV-GIR functions and long-run multipliers as will be explained in 
the following part. Last part of this section shows how to make statistical inference on 
time- varying linkage of stock markets. 

2.1 Stock Market Linkages and VAR Models 

Here we explain why we adopt a simple VAR model to analyze the stock market linkage. 
Let y^ denotes A;- dimensional vector representing the rate of return of k stocks at t. 
Preceding studies use VAR models to analyze linkage of stock markets, say VAR(p) model, 
for convenience of explanation: 

y^ = iy + Aiyt_^ + ■■■ + Apy^_p + Ut, t = 1, 2, ■ ■ • , (1) 

where i/ is a intercept term and Ut disturbance term at t, which is supposed to be inde- 
pendently identically distributed. 

Given the information set Xi_i available at t— 1, adapted to the equation ([1]), the follow- 
ing equation holds since I^-i is the cr-field generated from yi_i, 2/t_2, • • • and Ut-i, Ut-2 ■ ■ ■ ■ 

Vt - E[yt I It-i] =ut, t = 1, 2, ■ ■ ■ , 

since The equation above asserts that the excess returns of stocks at a certain period 
is "determined by accident" if the time series of stock returns follows such a stochastic 
process. This interpretation is closely related to t he efficient ma rket hypothesis (EMH) in 



the semi-strong form as well as predictability (see iFamal (Il970l )). We are quite interested 
in the time series of the excess returns of stocks. 

A convenient way to analyze it is expressing it as a moving average of past innovations, 
Wi's. Under a certain invertible conditon , we can convert the equation ([1]) to vector 
moving average form: 

yt = ^J' + ^oUt + ^lUt-l + $21*4-2, t = 1, 2, ■ • • , 

where $0 = ^k- When regarding the sequence {$j}^g as a function from t to k x k 
matrix $t, one usually call the secquence the impulse response function. Once the VAR 



coefficients are estimated from data, we can make some inferences such as the impulse 
response analysis and the Granger causality test to analyze the time series of the excess 
returns of stocks. 

For our purpose, analysis of time-varying structure of linkage of stock markets, we 
provide some new procedures and methods: a parameter constancy test for VAR model, 
estimation of TV- VAR model, computing TV-GIR functions, measuring long-run multi- 
pliers of TV- VAR. In the following subsections, we show the new procedures and methods 
making links to traditional procedures. 

2.2 Parameter Constancy Test for VAR Model 



There a re some param eter constancy tests, for example lAndrewsl (119931 ) and iNvblom 
(119891 ). JHanseni (119921 ) develops a parameter constancy test for linear and non linear 
models. The test is made under the null and alternative hypotheses: (Hq) the parame- 
ters are constant over time and (Hi) they follow a martingale process. In practice, we 
reformulate the equation ([T]) to extend Hansen's test for a VAR(p) model. 
Define the {k x T) data matrix as 



Y 



1 1/1,1/2, 



• Vt), 



where T is the number of time series observations. To reformulate the VAR(p) model with 
a intercept term into a regression formula, we introduce the following auxiliary vector and 
matrix: 

/ 1 \ 
Vt 



V Vt-p J 



and 



Z '■— {Zq, Zi,- ■ ■ , Zt-i) ■ 

Note that Zq is regaraded as a prior of estimation. Defining the {k x kT) disturbance 
matrix as 

we obtain a matrix form of the equation ([1]): 

Y = BZ + U, 

where a kx [1 + kp) matrix B is [u,Ai,- ■ ■ , Ap\. Using the vec operator, which transforms 
a (/c X T) matrix into an [kT) vector by stacking the columns, we have the following 
regression form of the equation ([T]): 



vec{Y) = {Z' ® Ik)vec{B) + vec{U), 



(2) 



where ® is the Kronecker product and Ik denotes k identity matrix. In case k = 2 and 
p = I, the equation ([2]) is specifically exhibibited as 
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In order to extend the procedure of iHansenl (119921 ). we rewrite the equation ([2]) as 

y = X/3 + n, 

where y = vec{Y), X = Z'^Ik and u = vec{U). Given y, we obtain the OLS estimate /3 = 
{X'X)~^X'y and the estimate of the covariance matrix E„ as UU' /T or UU' /{T — kp— 1) 
(see Liitkepohl (2005, ch.3) for details). 

Furthermore, to extend Hansen's (1992) procedure to our system, let us organize the 
above equation as follows: 



y, = Xt(3 + Ut, t = l,2,---,T. 



(3) 



For simplicity, we consider an OLS estimation of the equation (|3]) with the i.i.d. assump- 
tion, E[utUf] = T^uit = 1,2, ■■■. The OLS estimates (3 as well as Eu should hold the 
following equations. 



= J2Ker 



r=l 



and 



= y^{vech{ere'^) — vech{Tju)), 



T=l 



where et = y^ — Xtf3 and the vech operator is closely related to vec, which only stacks the 
elements on and below the main diagonal of a symmetric matrix. To construct Hansen's 
test statistic for a simultaneous linear equation system, define (fc(l + kp) + k{k + l)/2) 
vector, 

vec{X'^er) 
vech^ers'^) — vech(T,u) 



fr0X) 



r = l,2,--- ,T. 



Let fir denotes i-th component of /t(/3, S„), Sg = J2l=i fr and Sis = Sr=i fir ^^^ ^ ~ 
1, 2, ■ ■ ■ , k{l + kp) + k{k + l)/2, s = 1, 2, ■ ■ ■ , T. Note that we abbreviate f3, S„ when 
there is no confusion. Hansen's (1992) individual statistics of parameter constancy are 



Li 



1 ^ 
— V92 

TVi^ ''' 

T=l 



1,2,--- ,{k{l + kp) + k{k + l)/2), 



where Vi = ^^^^ /^v- His joint statistic is defined: 



1 ^ 



T = l 



where V = ^^^^ frfr- These test stati stics follow sin gular distributions represented by 
the Brownian motions and bridges. See iHansenl (1990) for the tables of the statistics. 



2.3 Non-Bayesian Time- Varying VAR Model 

If the parameter constancy test exhibits instability of VAR coefficients over time, we 
choose an alternative model that holds the alternative hypothesis on the parameter dy- 
namics and then estimate the model. Focusing our attention on linkage of stock markets 
which is supposed to vary over time, we suppose that only VAR coefficients vary over time 
while any components of the intercept term are invariant. In what follows, we present 
how to estimate the non-Bayesian TV- VAR model. Our approach has severa l good points 
in comparison with the Bayesian one (see, for example, ICogley and Sargentl (1200 ll . |2005| ) 
and IPrimiceril (20051)). 

We start by setting VAR coefficients varying over time for an ordinary VAR(p) model 
as follows: 

y* = ^ + Ai^tVt-i + ■■■ + Ap,tyt-p + Ut, t = 1, 2, ■ ■ ■ , T. 

This equation can be represented as follows: 



Vt 



u + AtZt.i + +ut, t = l,2. 



T 



(4) 



where 



At = [Ai^t ■ ■ ■ Ap^t] and Zt-i 



Vt-i 



Vt-p 



The equation (jl]) corresponds to an observation equation when we regard our TV- VAR 
model as a state space model. As to the corresponding state equation, let us assume the 
following random walk process: 



A 



i,t 



Ai., + Viu i = l,2, 



,p and t = 1, 2, 



T 



where each Vi^t is a fc x fc matrix of random variables, say, following normal distributions. 
Let Vt denote = [Vi^t ■ ■ ■ Vp,t]] Or equivalently 



vec{At) = vec{At^i) + Vt, t = 1, 2, ■ ■ ■ , T, 



(5) 



where each Vt is a k'^p-vector of random variables, veciVt) = vec{[Vi^t ■ ■ ■ ^,t]])- 

In place of the Kalman smoothing, we estimate v,Ai,--- , At by considering together 
the equations (jlj) and ([5]) as a simultaneous system of linear equations. The equation (jl]) 



turns to be 
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This is equivalent to the following equation. 



y = \l®h diag(Z^ ® Ik, Z( ® 4, ■ ■ ■ , Z^_i ® Ik, )] 

vec{[i^, Ai, A2,--- , At-i]) + vec{U), (6) 

where 1 = (1 1 ■ ■ ■ I)' G R . On the other hand, the state equation can be represented 

as: 
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where fec(Ao) is a prior vector of VAR coefficients, of whichever level we set do not 
influence our results except starting burn-in periods. 

Considering the following linear system, we estimate the coefficients of our time- 
varying VAR by using OLS. 
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where D = [l 4 | diag(Z^ ® h, Z[®Ik,--- , ^^-i ® h 
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The OLS estimate is: 
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As llto and Nodal (120121 ) show, the OLS estimate, /3, of the above hnear model provides the 
same estimation. At, (t = 1, 2, ■ ■ ■ , T) as those of the Kalman smoothing for the original 
state space model. In contrast to the traditional method of the Kalman smoothing, 
we can use familiar techniques of econometrics such as heteroskedastic autoeregressive 
consistent (HAG) estimations (see, for example. iNewey and WestI fll987l . Il994l )). Such a 
HAG estimate enables us to obtain the covariance estimate of E^,ec(At)5 (^ = 1; 2, ■ ■ ■ ,T) 
from the covariance estimate, S^. It is possible to construct the confidence intervals of 
the VAR coefficients for period by period by using the squared diagonal components of 
E^. Furthermore, as shown lA.ll in Appendix, our approach never confronts with the 
problem of identifiability unlike usual econometric analyses using least square techniques. 
In comparison with Bayesian approaches, our method is so simple that we are free from 
sometimes agonizing choice of prior distri bution of parameters ( for a t ypica l exar nple of the 
Bayse ian TV- VAR model approach, see 

fl2nn5h ). 



Gogley and SargentI (l200ll . |2005| ) and iPrimiceri 



2.4 Granger Causality for Non-Bayesian TV- VAR Model 

Similar to the impulse response functions, we have the Wald test statistic of Granger 
causality at t under the suppostion of locally stationary process. That is, we consider 
that it is possible to obtain the statistic at each period based on information available on 
the estimation of our TV- VAR model. 

The Granger causality test is essentially a test for a constraint of VAR coefficeients. 
Inspired by Liitkepohl (2005, ch.3.6), we generally consider the following tests for our 
TV- VAR model: 



Ho^t : C(Xt = c against iJo,i : Co^t 7^ c, (t = 1, 2, 



• T) 



where C is an (A^ x [k"^ + k)) matrix of rank A^, ctt is a parameter vector of TV- VAR 

coefficients at time t including drift terms or not, and c is an (A^ x 1) vector. Suppose 

that 

ctt- ctt \ d j^ I' f O \ ( Ea, O 

& - (T \\0 



-'at 

O E, 



holds. Then, 



\fT{Cctt - Cctt) 4 M (0, CE^/C) , t = 1, 2, ■ ■ ■ 
is derived. Hence, 

T{Cdtt - c){CT.^,C')-\Cctt - c)' A x\N). 
Resulting test statistic is: 

\wt = TiCat - c)iCJ:s.tC')-\Cctt - c)', t = 1, 2, ■ 
Using this statistic, we make a x^ test for the Granger causahty. 



.T. 
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2.5 Generalized Impulse Responses for Non-Bayesian TV-VAR 
Model 

In the previous section, we exhibit how to obtain the non-Bayesian TV-VAR estimate, 
At, {t = 1,2,--- ,T). We analyze time- varying structure of international stock market 
linkage based on the series of TV-VAR(l) coefficients. Summarizing his idea, he regards 
non-stationary time series generated by a state space model as patchwork made with 
locally stationary VAR models (see Kitagawa (2010, ch.l3) for detail). 

Here we show how to compute impulse response functions and related series of variables 
at each period t to elucidate the time- varying structure of linkage to stock markets. First, 
we give brief reviw of computing impulse responses of our TV-VAR model based on At. 
Given coefficients At of a VAR model, we can convert it to a vector moving average 
(VMA) representation under the condition that for any z such that |2;| < 1 

det(4 - Ai^tz Ap^tz^ ^ 0. 

Our idea relies on the generalized impulse responses, which are developed in order 



to a void the shorcoming whe n using the Cholesky decomposition (See lKoop et al.l ( 119961 ) 



and lPesaran and ShinI ( 19981 )). The generalized impulse response function Gly at horizon 



n of an underlying time series {Vt} is defined as: 

GIy{n, d, fit_i) = E (?;,+„ \ut = S, f]i_i) - E {y^^^ \ fi^.i) , (7) 

where Qt-i is the available information set at t. In our case our concern is the influence of 
a speciflc shock in a state equation (jlj) to state variables. We suppose that the two kinds 
of noises in (jl]) and ([5]) follow an i.i.d. distribution. That is, their covariance structure is 
represented: 

S„ O 

o s„ 

where S„ and S^, are covariance matrices of ut and Vt respectively. Furthermore, for 
convenience of calculation, we assume the noises follow the normal distribution with mean 
and the above covariance matrix. Under these assumptions, we can evaluate the two 
terms of right-hand side in ([7]). The reader should note that the difference between them 
only depends on shocks of the state equation (jlj) at t and that the shocks never cause 
variation of VAR coefficients At-^n for any n > since there is no correlation between Ut 
and Vt in equation (^ and ([5]) (see lA. 21 in Appendix for detail). This result alleviates our 
task to calculate the generalized impulse response function. 

Applying the above result and Liitkepohl (2005, ch.3) to our VAR(p) model allows us 
to have s'th step of the ordinary impulse responses, conditional on t— 1, $s j corresponding 
to At as follows: 

s 
$0,4 = h, ^s,t = Y^ ^s-j,tAj,t, s = 1, 1, 2, ■ ■ ■ . 

J=l 

Note that in case of p = 1 the above conversion rule becomes 

^s,t = Al s = 0,l,2,---. 



When we mention an impulse response function at t, $<, it means a mapping from non 
negative integer m to $m,t- 

Note that /x^ = E[yf.\Qt-i] varys over time in our framework in contrast to ordinary 
VAR models since 



fi = {Ik-A 



i,t 



A. 



p,t) 



\-i 



u. 



^it = u + Atu + Ail' + ■ ■ ■ , t = 1, 2, ■ ■ ■ , T. 

Employing the ordinary impulse responses at t and the following relation is essential to 
derive our TV-GIRs. 

We have 

(oo \ / oo 

^ ^i^tUt+n-z \Ut = d, Qt-1 ] - ^ [^ ^i,tUt+n-i \ ^t-1 
i=0 J \i=0 

= ($0,tWt+„ + ^l,tUt+n-l H 1- ^n,tS + ^n+l,tUt-l H ) 

- {^o,tUt+n + ^l,tUt+n-l H \- ^n,tUt + $„+l,iWt_l H ) 

Following the setting of shock by IPesaran and ShinI ( 19981 ). we obtain the unsealed gen- 
eralized impulse response funcion of the shock E{£t \ Ujt = 6j). Note that such a shock, 
which is also employed in the above derivation, is represented as 



Eie.\u,,^6,)^^l- £)(o)< 



jj^i 






a. 



u,]] 



0'u,Kj 





V J 



0"« 



jj 



where au,ij is the ij component of S,, . 

Following iPesaran and Shin! (jl998[ ). we choose the stength of shock: 



"7 ^u,ii- 



Then, the generalized impulse response of the effect of a shock in the j— th state equation 
at time t on y^.n i^ 



^ujj^n,t^uej. 



(9) 
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Finally, using diag(-), a diagonal operator whose value is the diagonal matrix which 
has the diagonal elements of the argument matrix or the argument vector on its main 
diagonal, we can derive the following representation of the generalized impulse response 
function. 



where T = S„ ■ diag(o-„ i\ a ^^ ^. ' 



u,KKJ 



2.6 Time- Varying Degrees of Market Efficiency and Propaga- 
tion of Market Fluctuations 

Thus, the excess returns of stocks at t, Zt = tft — A^i, is of interest to us when we analyze 
the time-varying structure of stock markets. In case of the efficient market hypothesis, 
we are interested in whether E[zt \It-i] = 0,t = 1,2,- ■ ■ are independently distributed 
or not. We could analyze it by examining the impulse responses at t, say, whether the 
s'th responses to an impulse are zero or not for any s > 1. 

We regard such time dependent variables of a TV-VAR model as of a locally stationary 
model. These impulse responses would give us information on how a shock appeared 
in a market at a period propagates over time locally. Furthermore, we are concerned 
about total effects of such an exogenous shock to each variable, which is calculated by 
cummulated summation of ^s,t if shocks are uncorrelated to each other. Specifically, 

s 

Taking limit of s to oo, we have 

oo 

i=i 
If Jfc — Ai^t — • • ■ — Ap^t is nonsigular, the following equation holds for VAR(p). 



^oo,t = (4 - ^i,t Ap^t 



,-1 



Then, we use this relation to obtain total effects, '^s,t for each t. We call '^s,t and ^oo,t 
the long-run multipliers, which negatively reflect the degree of market efficiency in the 
semi-strong sense. That is, the larger the value of long-run multipliers are, the more 
inefficent the stock markets are. 

In order to examine the propagation pattern of a certain market shock to both the 
market itself and other ones, we consider the underlying shock that is uncorrelated to any 
shocks in other markets. In other words, responses to an indpendent impulse are signifi- 
cant. When a VAR model is employed, however, the estimated covariance matrix of the 
error term, S„, is not always diagonal. Conventionally, the use of Cholesky decomposition 
has been suggested (see, for example, Liitkepohl (2005, ch.3)). However, this approach is 
not invariant to the ordering of variables in a VAR models. We employ the generalized 
impulse responses for a VAR model to avoid this shortcoming. 
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Liitkepohll ( 120051 ) shows that an orthogonahzed impulse response function of a VMA 
model is obtained simply by multiplying each coefficient matrix by P, the Cholesky de- 
composition matrix of the covariance matrix of shocks. In place of P, we use T defined 
in the previous subsection to obtain the generalized impulse response function. Thus, we 
define their corespondences to $5 j, \E's j and \E'oo,t for the generalized impulses responses 
as follows. 

e,,t = <l>,,iT, Es,t = ^s,tT, Soo,i = ^oo,tT. 

Practically, we estimate for the matrix T by using T to obtain the estimates, 6^,4, S^ j 

and Soo,t. 

We exhibit the relations among the variables: 




The reader should notice the significance of each estimates in the above diagram. In 
particular, the TV-GIRs Qt and the time-varying long-run multipliers matrices Sqo t play 
a very vital role in our paper. The former provide us with information about how a market 
shock infiuences on the other market over time. It concerns the relation between two 
markets through an instaneous shock in a market. The latter teaches us the time- varying 
structure of the relation between two markets from the point of view of information 
efficiency. Specifically the diagonal compon ents of the matrix at a time Hoo,t exhibit 
market inefficiency in the semi-strong sense of lFamal ( 1l970l . Il991l ). The non-diagonal ones 
refiect how strong a market propagates its information to another market. We call the 
non-diagonal components of Hoo,t at time t the degree of propagation]^ 



3 Data 

This paper uses the monthly returns for the S&P500 and Tokyo stock price index (S&PSOO 
and TOPIX) from January 1955 to September 2012, obtained from Robert Shiller's web- 
site (S&P500) and The monthly statistics report of Tokyo Stock Exchange (TOPIX). In 
practice, we take log first difference of time series of the stock price index. Figure [T] 
presents time series plots for each variables. 

(Figure [T] around here) 

For the estimations, all variables that appear in the moment conditions should be 
stationary. To check whether the variables satisfy the stationarity condition. We use the 



^Following the method of Liitkepohl (2005, ch.3, Proposition 3.6), we can compute the confidence 
bands of the TV-GIRs and their related variables (see lA.2l in Appendix for detail). 
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ADF-GLS test of Elliott et al.l ( 1l996l ). Table [T] provides some descriptive statistics and 
the results of the ADF-GLS tests. For all the variables, the ADF-GLS test rejects the 
null hypothesis that the variable contains a unit root at conventional significance levelso 

(Table [1] around here) 

4 Empirical Results 

First, we present a preliminary VAR estimation with whole sample and examine whether 



the VA R estimates change over time or not using the paramter constancy test of [Hansen 



(Il992l ). Second, we estimate the our TV- VAR model, and calculate Granger causality test 
statistics as shown in section |2j Third, we also estimate the long-run multipliers which 
are based on the TV-GIRs to measure the degree of mark et effic i ency and pr opagation of 



market fluctuations in the sense of semi-strong form (see, iFamal (jl970l . Il99ll ) for details). 



4.1 Preliminary Estimation 

We assu r ne a m odel with constant and use the Bayesian information criterion (BIG) of 



Schwarzl (119781 ) as an optimal lag order selection criteria in a VAR estimation. In our 
estimation, we choose a first-order vector autoregressive model. Table [2] shows that the 
preliminary results for an VAR(l) model using the whole samplejj 

(Table [2] around here) 

The VAR estimates are statistically significant at conventional levels except for first- 
order vector autoregressive term, R'^^^ , in the equation R^^ . And the results of Granger 
causality test show that the returns on U.S. stock market infiuence the returns on Japanese 
stock market during the above sample period, but not vice versau 

Now, we investigate whether there exhibits a parameter constancy in the above VAR 
model using Hansen's (1992) tests under random parameters hypothesis. Table [2] also 
shows the results of Hansen's (1992) joint parameter constancy test with variance, the 
VAR estimates are simultaneously time- varying, we reject the null of constant parameters 
against the parameter variation as random walk at the 1% statistical significant level. 
Therefore, we estimate the time- varying parameters of the above VAR model to investigate 
whether gradual changes occur in the U.S. and the Japanese stock market. 



^We confirm that there are no size distortions that Elliott et al.l (jl996h and iNg and PerronI ( 2001r ) 



pointed out in making the ADF-GLS test for small samples (see a column of the table [T] for details). 
Therefore, we use the Modified Bayesian Information Criterion (MBIC), not the Modified Akaike Infor- 
mation Criteria (MAIC) to select an optimal lag order for the ADF-GLS tests. 

•^We use the heteroskedast icity and autocorrelation consistent (HAG) covariance matrix estimator of 
Newev and \\tetl (Il987l.ll994l). 



iToda and Yamamotd ( 1995 ) demonstrate that the Granger causality test should be carefully made 
for VAR estimation formulated in levels, since the test concerns general restrictions on the parameter 
matrices even when the processes may be integrated or cointegrated of an arbitrary order. Their result 
has nothing to do with our approach because we do not use level variables and we choose one as the lag 
order of VAR using BIG. 
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4.2 Non-Bayesian TV-VAR Estimation and Granger Causality 
Test 

We adopt new estimation methods as shown in section 2.3 for our TV-VAR(l) modeL 
Figure [2] exhibits the time- varying estimates of the VAR(l) modeL 

(Figure [2] around here) 

We confirm that all first-order vector autoregressive terms change over time in our TV- 
VAR(l) model. The esti mates show that the struct ure of stock markets fluct uates wildly 



through time as shown in llto and Sugiyamal (120091 ) and llto and Nodal (120121 ). 



Figure [3] shows that the results of the Granger causality test which are based on the 
TV-VAR(l) estimates. 

(Figure [3] around here) 

We confirm that the results of Granger causality test change over time. And those results 
are quite different from when we employ a time-invariant VAR(l) model. In particular, 
we find that: (1) the top panel in figure [3] provides the Granger causality from the returns 
on U.S. stock market to the returns on the Japanese stock market in two decades from 
1970 and in several years from late 1990s, and (2) the bottom panel in figure [3] displays 
the Grenger causality from the returns on Japanese stock market to the returns on the 
U.S. stock market in the early 1980s and the early 1990s. 

4.3 Time- Varying Generalized Impulse Responses and Long- 
Run Multipliers 

In this subsection, we estimate the long-run multipliers which are based on TV-GIRs to 
measure the degree of market efficiency and propagation of market fluctuations in the 
sense of semi-strong form. And we analyze the total effects of exogenous shocks. 

Figure H] shows the three-dimensional representation of the TV-GIR functions in our 
TV-VAR(l) model. 

(Figure H] around here) 

We find that the time paths of effects of an exogenous shocks to the each stock mar- 
kets vary widely with time. And we estimate the long-run multipliers which are based 
on the TV-GIRs to measure the degree of market efficiency and propagation of market 
fluctuations. 

The diagonal panels in figure [5] present the degree of semi-strong form of market 
efficiency. 

(Figure \5\ around here) 

We confirm that the degree of semi-strong form of market efficiency changes over time 
for each stock markets. The degree of market efficiency declines (rises) in the period 
of bubbles (panics: market crashes and economic crises) in both the stock markets as 
is obvious from the table [31 In periods of bubble with full of unreliable information, 
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some investors could gain huge windfalls, which suggests that the market is inefficient; in 
periods of panic, an average investor hardly beats the stock market overall, which suggests 
that the market is efficient. 

(Table [3] around here) 

And the degree of market efficiency in the semi-strong sense shows that the U.S. stock 
market is more efficient than the Japanese one at most periods. 

The non-diagonal panels in figure [5] exhibit absolutely different degree of propagation 
of market fluctuations in U.S. and Japan. According to the deviations from zero in the 
two panels, we find several facts. First, the degree of propagation of market fluctuations 
in U.S. is two times larger than Japanese ones at most periods. In other words, the U.S. 
stock market is hardly influenced by the persistent shocks from the Japanese stock market, 
but the Japanese stock market recieves a great amount of influence from the U.S. stock 
market. Second, the U.S. stock market strongly influenced the Japanese stock market 
in the periods of "1973 Oil Crisis" and "Booming 1980s". Third, the events in the mid 
2000s, "U.S. Housing" and "Bankruptcy of Lehman Brothers", have smaller influence on 
the Japanese stock market than in the above two eras. 

5 Concluding Remarks 

In this paper, we develop new methodologies to analyze the time-varying structure of 
international linkages and market efficiency in stock markets. In particular, we ffist es- 
tablish whether the parameters of the standard VAR model are time- invariant or not. The 
results of Hansen's (1992) test indicates possibility that those parameters are not time- 
invariant. Then, we estimate our non-Bayesian TV- VAR model and calculate Granger 
causality test statistics. And we estimate the long-run multipliers which are based on 
the TV-GIRs to analyze the total effects of exogenous shocks as the degree of market 
efficiency and propagation of market fluctuations in the semi-strong sense. Our empirical 
results provide a new perspective that the stock market linkages and market efficiency in 
the semi-strong sense change over time. 

First, we find that the international stock market linkages between the two countries 
vary with time from the following two evidences. The results of the time- varying Granger 
causality tests show that the strength of the Granger causahty depends on period by 
period. Moreover, the degree of propagation of market fluctuations in U.S. is two times 
larger than Japanese ones at most periods. It indicates that the international stock market 
linkages between the two countries changes considerably for some periods, for instance, 
the Japanese stock market is not so much greatly influenced by the bankruptcy of Lehman 
Brothers as the U.S. ones. Our findings show the international linkage of stock market 
varying. While the preceding studies certainly approve the existence of the linkage, we 
stress strength of the linkages differ periods by periods. 

Second, we also find that the degree of market efficiency in the semi-strong sense 
changes over time for each stock markets. The degree of market efficiency declines (rises) 
in the period of bub bles (panics: rnarket crashes and economic crises) in both the stock 



markets as same as llto and Nodal (120121 ) pointed out. In periods of bubble with full of 
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unreliable information, some investors could gain huge windfalls, which suggests that the 
market is inefficient; in periods of panic, an average investor hardly beats the stock market 
overall, which suggests that the market is efficient. And the degree of market efficiency in 
the semi-strong sense shows that the U.S. stock market is more efficient than the Japanese 
one. 

At last, as pointed out in the last part of the previous section, our empirical results 
well correspond to the major historical events in international finance in this half century. 
Although varying stock market linkage certainly reflects increased structural changes in 
our economy, what else shock affects the stock markets is still an open question. Each 
country's economic policy, say, its flscal policy, would affect the international flnancial 
market. This issue is left for our future research. 



Appendix 



A.l Regressor of State Space Regression Models 

This section shows the reader a signiflicant property of the method that we present for 
dealing with a TV-VAR(l) model in section |21 The method enables us to handle the 
model as a conventional econometric model by identifying the corresponding state space 
model with a linar regression model. Recall that the observation equation and the state 
one appeared in the state space model in section [2j the equations, (jl]) and ([5]). 

We prove that estimating method of the smoother of a state space model guarantees 
the exact estimate whatever data is applied. Mathematically we assert that the regressor 
matrix of our underlying linear system has full rank. This implies that our estimation 
ne ver suffers fr o m mu lti colinearity and that our model is always identiflable in the sense 
of iRothenbergI fll97ll ). The following results do not depend on whatever time- varying 
model is assumed; they cover our TV-VAR model, a typical state space model. 

The model consists of two equations as follows: 



Vt = ^tPt + 'i^u 



«/~Ar(0,S, 



t 



1,2, 



T 



(A.l) 



and 



l3,^-, = Bt(3, + Vt, i;/~Ar(0,E,), t = l,2,---,T. (A.2) 

Note that 7//s are /c— dimensional vectors and /^^'s are m— dimensional vectors as para- 
maeters. Xt and Bt are k x m and m x m matrices respectively. Notice that two dis- 
turbances, Ut and Vt, follow iid noromal distributions with zero means and covariance 



matrices, 5]^ and T,^, respectively. 



By transformig the equation (I A. 2 p to the following form, we can regard the state equa- 
tion as a linear regression model. 







iid 



Btf3t-l(3t^, + vt, t;/~Ar(0,S„), t = l,2,---,T. (A.3) 

Combining (lA.ip with (IA.3P into the following linear system, we can employ conventional 
econometric techniques to obtain the Kalman smoother. 



y 



y 

7 



X 

w 



/3 + 
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where 
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Our main theorem of this appendix is as follows: 

Theorem 1 The regressor matrix has full rank. That is 



B' 



T-l 



rank 



X 
W 



mT. 



Proof: Let X denote the regressor matrix. It is of size {k + m)T x mT. Thus, 
rankX < mT. Clearly, each ranks of the two submatrices, W and X, are fewer than or 
equal to mT. Formally, 

rankX < mT, 



and 



rankVF < mT. 



On the other hand, rankVF = rankVF'. Since W has the following reduced row 

echelon form: 

/-IB, O \ 



-I 



\ o 



Bt~i 
-/ 



its rank is clearly mT. Finally, mT < rankX < mT. This implies rankX = mT. D 

The reader should notice that the rank of X depends neither the data matrices, 
Xi, X2, ■ ■ ■ , Xt, nor the state transition matrices, Bi, B2,- ■ ■ , -Bt-i- This signifies that 
classical least square techniques such as OLS or GLS are always applicable to estimate 
linear time-varying models that can be represented as state space models. 

In this paper, our TV-VAR(l) model with constant drift. Such a model has two 
types of parameters to be estimated: time- varying and time- invariant. In the rest of this 
appendix, we confirm that the above discussion holds even if linear time-varying model 
such as our model has time-invariant parameters. In the case, we should modify the 
equation (]A.1|) as follows: 



yt 



ZtOL + Xi/3j + Ui, 



w/~Ar(0,S, 



t = l,2, 
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where a is a £— vector of time invariant parameters and Zt is a k x i matrix of data. The 
regressor matrix turns to be 

O W 

where the above O is mT x i matrix of zero and 

Z2 

\Zt J 
The following theorem holds. 

Theorem 2 

rank 



Z X 

o~w 



lankZ + rriT. 



Proof: Let X denote the regressor matrix in the case. According to Theorem [H rankX 
is mT. Thus, using the reduced row echelon form of W, say, Wq, we can transform X 
to the following form: 

O 

Wo' 

Then there are two non-singular matrices, P of size {{k + m)T) x ((/c + m)T) and Q of 
size {£ + mT) x {£ + mT), such that 



PXQ 



Considering rankZ = mT and 



rank 



Z 

o 



z o 

O Wo 



rankZ', 



we have rankX = TenakZ + m,T. D 

Note that we can deal with the case of the TV-VAR(l) with time-invariant drift term 
such as our model in this paper. As appeared in section 12. 3[ one regards Z as 

//\ 
/ 

where each / denotes kx k identity matrix. Then, rankZ = k. Considering the number of 
parameters to be estimated is fc + fc^T and m, = k"^, our TV-VAR model has no problem 
of identification. 
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A. 2 Generalized Impulse Responses and Their Confidence Bands 

This section presents several mathematical propositions on which our empirical analyses 
rely. First, we address the important property of the generalized impulse responses for 
our TV-VAR model. That is, they are derived from the VMA model such that we set 
shocks concerning the time-varying parameters of our TV-VAR model zero for periods 
foUwoing the starting period. Then, we employ the matrix T in section 1231 in place of the 
orthogonal shocks derived from the Cholesky decomposition. Second, we can compute 
the confidence bands of the genenralized impulse responses and other our estimates since 
they depend only on the estimates of the TV-VAR coefficients and covariance matrix of 
the disturbance terms of the model. 

Here we show the first proposition with respect to the generalized impulse responses. 

Proposition 1 Consider TV-VAR(p) model, the equation ^ and ^. Suppose that 



(Ut,Vt 



-ATjO, 



-'u 

O 



o 



The generalized impulse response function of y^ for TV-VAR (p) at horizon n is obtained 
from a VMA(oo) model conditional on Vtt-i derived by setting = Vt = Vt^i = Vt^2 ■ ■ ■ 
for the TV-VAR(p) model. 

Proof: Fix t and recall that 



Vt+i = At+iZt-i+i + Ut+i i = 
vec{At+i) = vec{At^i+i) + Vt+i i 



1,2,---, 



1,2,---, 



and 



Vt+i 



iid , r 








s„ o 

o s„ 



Thus, following iKoop et al.l (119961 ). we have the shocks to obtain the generalized impulse 
responses as follows: 

E{et I £jt = i^uj, 0) , fii-l) = {o-ulj, <7u2j, ■ ■ ■ , O'uKj, 0, ■ ■ ■ ,0) (Tujj^uj 



s„ o 

o s„ 



^u^ujO'ujj 



-U] 



o 



o 

6 



a, 



-^6 ■ 



UJI 



where auij is the ij component of H^- This fact allows us to ignore the shocks concern- 
ing time-varying parameters in our TV-VAR model when we derive the corresponding 
VMA(cx)) model. D 

The above proposition shows that the generalized impulse responses can be obtained 
by ignoring shocks on time-varying parameters. Thus, we can apply similar derivation 
process of impupulse responses from orthogonalized shocks in Liitkepohl (2005, ch.3). We 
emplooy the matrix T in section [2751 to obtain the generalized impulse responses and other 



19 



stats derived from them as pointed in section [2j The remainder of this section shows how 
to compute the confidence bands of the TV-GIR functions and their related variables. We 
can construct the confidnence band under the assumption of asymptotic normality of the 
parameters estimated, vec{Ai, A2, ■ ■ ■ , A^) for TV-VAR(p). Since all we need are not the 
(k'^pT) X (k'^pT) covariance matrix of vec{Ai, A2, ■ ■ ■ , A^) but the T covariance matrices 
of vec{At), (t = 1, 2, ■ ■ ■ , T) as well as the covariance matrix of cr = vech{T,u)- We take 
T submatrices with size k^p x k'^p from the diagonal positions of vec{Ai, /I2, ■ ■ ■ , At) as 
the covariance matrices of vec{At), t = 1, 2, ■ • ■ , T. Let Sq,j denote the covariance matrix 
of vec{At), t = 1, 2, ■ ■ ■ , T in what follows. 

With respect to the covariance matrix of cr = fec/i(Sij), noted So-, the following 
formula holds. 



2£)+(S, 



^u)iD 



+\' 



where D^ = [D'^^Dk) ^D'k^ ^^at is the generalized inverse of dupulication matrix Dk, 
which is a constant matrix. In case of fc = 2 and p = 1, 



D, 



/ 1 \ 
1 
1 

Vo 1/ 



£> 



10 

i I 

1 



We have its estimate, Sa- setting S^^ = UU' /T, where U is the matrix of our OLS residuals 
(see Liitkepohl (2005, p.93) for detail). 

Here we present a theorem to summarize the asymptotic properties of our featured 
statistics under the supposition of Proposition 10 



Theorem 3 Suppose that 



oct - at 
& — cr 

d 



Am 



o 
o 



So, o 

o s^ 



where 



Vf{^s,t - $s,t) ^ AT (0, a,tE^,G; J , s = 1, 2, ■ ■ ■ , t = 1, 2, 



■ T 



G 



Ai,t A2,t 

h o 

o h 

o o 



- m,ty 



m=l 



Ap-l^t Ap^t 

o o 



h o 



Ut 



Ut 





kpxkp 




and {k x kp) matrix J = [Ik O ■ ■ ■ O]. 

^This part greatly owes to the results of Liitkepohl (2005, cli.3. sec. 7). 
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VT(^n,t - ^n,t) A M (0, F„,iS^,F;,) , s = 1, 2, ■ ■ ■ , t = 1, 2, ■ ■ ■ , T, 



where 



Fn,t — G 



i,t 



+ a 



n,t- 



// (/fc — At) zs nonsingular, 



VT(^oo,t - ^oo,t) A M (0, Foo,iS,,,F^ J , s = 1, 2, ■ ■ ■ , t = 1, 2, ■ ■ ■ , T, 



where 



i^oo,t = $'oot®'^oo,i- 



\/r(0,,, - 0,,i) A u (0, c,,iSa,c;, + c',,tS^c';,) , s = i, 2, ■ ■ ■ , t = i, 2, 

where 

Co = 0, Cs,t = {r'®h)Gs,t, s = l,2,---, 

Cs,t={h^'^s,t)H, s = 0,1,2,- ■■ , 



anc? 



da' 
( o'J 



da' 



B(a) 
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\ ^,7 / 

\/T(S,,i - S,,i) A AT (0, Bn,t^s.A,t + Bn,t^^B'^^,) , s = 1,2,- ■ ■ , t = 1,2, 



T 



where 
and 



5„,t = (T' ® 4)F„,i, n = l,2, 
Bn,t = {h®'^n,t)H, n = l,2. 



T 



VT{Eao,t — 2oo,t) — > A/" (O, Boo,t^a.tB'oo^t + -^oo,tSo-B^ J 



s = l,2,---, t = l,2,---,T, (A.4) 



where 
and 



B^,t = {r'®h)F^,t, t = l,2. 



Boo,t = {h 8) ^oo,t)i^, t = 1, 2, 
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Proof: With respect to ^s,t, ^s,t and \l/oo,t, our derivation directly follows the proof of 
Liitkepohl (2005, ch.3, Proposition 3.6) based on the delta method. Recall that 



T = Sii • diag n/A/diag(Si 
The key relation follows from the delta method. 

\/Tvec[Qi^t - Qi,t) ^ A/ 0, — S^, \ -— S^- ^ . 

\ ool' OCX. ocr' OCT J 



(9t;ec($i,iT0 = ^y' ^ j ^^vec{^i^t) 



da' 



i-l 

T/ /|'\i— 1— m ,0, <T> 



(9t;ec($i,iT) 9 



da- , dcr' 



\m=C 
/ \ 

S^,diag(l/diag(S^,)) 



"^/^ 
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\ B((T) J 



C<T Ccr 



a/T (^S„,t - En,t) 4 AT (0, 5„,iSA,5;^, + Bn,t^^B',J 
dveci^n,tr) 



B. 



n,t 



da 



dvec{'^n,t) 



da' 

da' 
,r.r,, ^ ^ v^ dvec(^it) 



= (r®4)5^a 
= (r 4)F„,i, 

where F„,i = Gi,i + G2,t H h Gn,t 

dveci^n^tT) 



Bn,t 



(4®^ 



dvec{T) 



n,t 



dcr' 
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Soo,t = (4 - Ai,t Ap,t)-^T{&), 

Foo,t = (^oo,t, ■ ■ ■ , ^oo,t) ®^oo,i, 

^ V ' 

P times 
yT f Soo,t — 2oo,t j — > A/" (O, Boo,t^atB'oo,t + Boo,t^&Boo,t) , 

where 

5oo,t = (T' ® 4)F^,i, 

and 

— OvccCY) 

Boo,t = (4 8) ^oo,t) ^ , ■ 

Note that ^s,t is independent of cr. D 
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Figure 1: The Returns on U.S. and Japanese Stock Market 
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Table 1: Descriptive Statistics and Unit Root Tests 





Mean 


SD 


Min 


Max 


ADF-GLS 


Lag 


<P 


M 


R^^ 


0.0054 


0.0360 


-0.2280 


0.1135 


-18.4965 





0.3366 


692 


^JPN 


0.0046 


0.0430 


-0.2439 


0.1336 


-18.8900 





0.3183 


692 


Notes: 



















(1) "/jf-^'S'" denotes the returns for the S&P500 stock index, and "/J-^-P^" denotes the retunrs for 
the TOPIX stock index. 

(2) "ADF-GLS" denotes the ADF-GLS test statistics, "Lag" denotes the lag order selected by the 
M BIC, and "(/>" denotes the coefficients vector in the GLS detrended series (see equation (6) 



Ng and PerronI ()2001l) l 



(3) In computing the ADF-GLS test, a model with a time trend and a constant is assumed. The 
critical value at the 1% significance level for the ADF-GLS test is "—3.42". 

(4) "A/"" denotes the number of observations. 

(5) R version 2.15.2 was used to compute the statistics. 
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Table 2: Preliminary Estimation 



Constant R[ 



UB' 



R-i 



7TW 

1 



GC 



K' 



'C 



i?r 



J^JPN 



0.0040 
[0.0015] 
0.0028 
[0.0016] 



0.2568 
[0.0419] 
0.0827 
[0.0476] 



-0.0200 

[0.0327] 
0.2932 
[0.0411] 



3.1823 

(0.0747) 

0.3757 

(0.5400) 



0.0586 



0.1052 



91.6098 



Notes: 



(1) "R^^i" denotes the lagged returns for the S&P500 stock index at time t — 1, and 
"-R/-?/^" denotes the lagged retunrs for the TOPIX stock index at time t — 1 . 

(2) "GC" denotes the Wald statistics for the Granger causality test. The null hy- 
pothesis that U.S. (Japan) do not Granger-cause Japan (U.S.). 

(3) "^2" denotes the adjusted R"^, and "Lc" denotes the Hansen's (1992) joint L 
statistic with variance. 

(4) Newey and West's robust standard errors and p- values for the Granger causality 
test statistics are in brackets and parentheses, respectively. 

(5) R version 2.15.2 was used to compute the estimates and the statistics. 
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Figure 2: TV-VAR Estimation 



US to US 



Japan to US 



O 
O 

cr 



oq 
o 



o 



o 
o 






in 


(0 


o 






c- 




(U 




o 




% 


o 


o 


o 


o 




rr 




<■ 




> 


IT) 




o 




US to Japan 



Japan to Japan 





tn 


to 


o 






c- 




(U 




o 




f, 


o 


o 


o 


o 




rr 




<r 




> 


IT) 




o 




O 
O 

cr 



oq 
o 



o 



o 
o 




Notes : 

(1) The dashed hnes represent upper and lower bounds of the time- varying 
VAR estimates at the 95% intervals for each period. 

(2) The Newey and West's HAC estimator is used for calculating the asymp- 
totic confidence intervals. 

(3) The dotted lines also represent the VAR estimates with whole sample. 

(4) R version 2.15.2 was used to compute the estimates. 
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Figure 3: Time- Varying Granger Causality Tests 
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The dashed hne denotes asymptotic critical velues at 10% significant levels 
for the Granger causality test statistics. 
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Figure 4: Time- Varying Generalized Impulse Response Functions 
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Long-Run Multipliers 



Long-Run Multipliers 
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Table 3: Dates for Market Crashes, Economic Crises, and Bubbles 



Events 


Start Dates 


End Dates 


Market Crashes 






Black Monday 


1987:10 


- 


Bankruptcy of Lehman Brothers 


2008:10 


- 


Economic Crises 






1973 Oil Crisis 


1973:10 


- 


1979 Oil Crisis 


1979:01 


- 


Asia Financial Crisis 


1997:07 


- 


Bubbles 






U.S. Booming 1980s 


1982:12 


1990:07 


Japanese Asset Price 


1986:12 


1991:02 


Dot-com 


1995:08 


2000:03 


U.S. Housing 


2001:04 


2007:10 


Notes: 







(1) The end dates for market crashes, and economic crises are unspecified 
because of no consensus. 

(2) The dates for bubbles in U.S. basically follow the NBER business 



cycle report (http://www.nber.0rg/cycles.h.t1nll. 



(3) The dates for bubbles in Japan basically follow the index of 
business condition by the Cabinet Office, Goverment of Japan 
(http: //www. esri . cao .go . jp/en/stat/di/di-e .html). 
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